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We perform a joint analysis of the data from PHENIX at RHIC and ALICE at LHC on the 
nuclear modification factor Raa- The computations are performed within the light-cone path 
integral approach to induced gluon emission. Our results show that slow variation of Raa 
from RHIC to LHC energies indicates that the QCD coupling constant is suppressed in the 
quark-gluon plasma produced at LHC. 



1. One of the most striking results of experiments at RHIC is strong suppression of high-p^^ 
hadrons in ^^-collisions ^ (called "jet quenching"). Recently, a similar effect has been observed 
in the ALICE experiment at for Pb + Pb collisions at ^/s = 2.76 TeV. The most natural 
reason for this phenomenon is parton energy loss (radiative and collisional) in the hot quark- 
gluon plasma (QGP) produced in the initial stage of AA-collisions. It is of great interest to 
perform a joint analysis of the RHIC and LHC data. It is interesting since variation of the 
nuclear modification factor Raa from RHIC to LHC energies should not be very sensitive to 
the systematic theoretical uncertainties that are rather large. These uncertainties come mostly 
from multiple induced gluon emission. The available theoretical approaches to radiative induced 
gluon emission 3 4 5 6 7 8] ^j.^ restricted to one gluon emission, and the multiple gluon emission 
is usually evaluated in the approximation of independent gluon radiation^Sl. 

In this talk, I will present results of an analysis of the data on Raa for Au + Au collisions 
at y/s = 200 GeV from PHENIxHSl 

and for Pb + Pb collisions at y/s = 2.76 TeV from ALICE 
121. The analysis is based on the light-cone path integral (LCPI) approach SI. We evaluate the 
nuclear modification factor using the method developed in^J^. A major purpose of this analysis 
is to decide whether the variation of Raa from RHIC to LHC indicates that the QCD coupling 
constant becomes smaller in the plasma produced at LHC, which is hotter than that at RHIC. 

2. The nuclear modification factor Raa for a given impact parameter b can be written as 



Here p-p is the hadron transverse momentum, y is rapidity (we consider the central region 
y = 0), TAAib) = J dpTA{p)TA{p — b), Ta is the nucleus profile function. The differential yield 
for high-pj- hadron production in AA-collision can be written in the form 
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where dam{N + N ^ h + X) / dprpdy is the medium- modified cross section for the N+N — )• h+X 
process. Similarly to the ordinary pQCD formula, we write it as 

dam{N + N ^h + X) dz da{N + N^i + X) 

d^ -^Jo ■ (3) 

Here = p^/z is the parton transverse momentum, da{N + N ^ i + X) / dpljdy is the hard 
cross section, -DJ^^ is the medium-modified fragmentation function (FF) for transition of a parton 
i into the observed hadron h. For the parton virtuality scale Q we take the parton transverse 
momentum pj,. We assume that hadronization occurs outside of the QGP. For jets with E < 100 
GeV the hadronization scale, /i/j, is relatively small. Indeed, one can easily show that the L 
dependence of the parton virtuality reads Q'^{L) ~ max (Q/L, Qq), where Qo ~ 1 ~ 2 GeV is 
some minimal nonperturbative scale. For RHIC and LHC, when tqqp ~ Ra (tqgp is the typical 
lifetime/size of the QGP, Ra is the nucleus radius), it gives fj-h ~ Qo (for E < 100 GeV). Then 
we can write 

dz' 

D^/,{z,Q)^ —D,/^{z/z',Qo)Dj)^{z',Qo,Q), (4) 

where Dfi/j(z,Qo) is the vacuum FF, and DY/^{z' ,Qo,Q) is the medium-modified FF for tran- 
sition of the initial parton i with virtuality Q to a parton j with virtuality Qq. For partons 
with E < 100 GeV the typical length scale dominating the energy loss in the DGLAP stage is 

relatively small ~ 0.3 — 1 This length is of the order of the formation time of the QGP 

To ~ 0.5 fm. Since the induced radiation stage occurs at larger length scale / ~ tq tqgp, to 
the first approximation one can ignore the overlap of the DGLAP and induced radiation stages 
at a Then we can write 
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T/,{z,Qo,Q)= / —Df/f{z/z',Ei)Df}f^^^^{z',Qo,Q), (5) 



where Ei = Qz', DJ^f is the induced radiation FF (it depends on the parton energy E, but not 
virtuality), and Df^jf^^^ is the vacuum DGLAP FF. 

We have computed the DGLAP FFs with the help of the PYTHIA event generator^IZ] One 
gluon induced emission has been computed within the LCPI formalism SI using the method 
elaborated in^l^. As in^ISIIl] -^g take mq = 300 and mg = 400 MeV for the quark and gluon 
quasiparticle masses. Our method of calculation of the in-medium FF via the one gluon prob- 
ability distribution is described in detail in and need not to be repeated here. We just 
enumerate its basic aspects. The multiple gluon emission is accounted for employing Landau's 
method as in^. For quarks the leakage of the probability to the unphysical region of /S.E > E 
is accounted for by renormalizing the FF. We also take into account the q ^ g FF. Its nor- 
malization is fixed from the momentum conservation for q ^ q and q ^ g transitions. The 
normalization of the g — )• q FF is also fixed from the momentum sum rule. The collisional 
energy loss, which is small is taken into account by renormalizing the temperature of the 
QGP for the radiative FFs using the condition: AEradiT^) = AEj.ad{To) + AEcoi{To), where 
^^rad/coi is the radiative/collisional energy loss, Tq is the real initial temperature of the QGP, 
and Tq is the renormalized temperature. 

We calculate the hard cross sections using the LO pQCD formula. To simulate the higher 
order i^T-factor we take for the virtuality scale in the value cQ with c = 0.265 as in the 
PYTHIA event generator We account for the nuclear modification of the parton densities 
(which leads to some small deviation of Raa from unity even without parton energy loss) with 
the help of the EKS98 correction!^. For the vacuum FFs we use the KKP parametrization^ISl, 

As in we evaluate the induced gluon emission and the collisional energy loss for the 
running frozen at some value q{'' at low momenta. For vacuum a reasonable choice is 
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ai' 0.7. This value was previously obtained by fitting the low-x proton structure function F2 



within the dipole BFKL equation^J^. To study the role of the in-medium suppression of Og we 
perform the computations for several smaller values of . 

3. We describe the QGP in the Bjorken model which gives T^To = T^T. We take tq = 0.5 
fm. To simplify numerical computations for each impact parameter b we neglect variation 
of the initial temperature Tq in the transverse directions. We evaluate its value using the 
entropy/multiplicity ratio dS / dy j dN^h/ drj ~ 7.67 obtained in^l^. For the central Au + Au 
collisions at ^ = 200 GeV Tq 300 MeV and for Pb + Pb collisions at ^/s = 2.76 TeV Tq ^ 400 
MeV. For the nuclear density we use the Woods-Saxon nucleus density with parameters as in 

The fast parton path length in the QGP, L, in the medium has been calculated according to 
the position of the hard reaction in the impact parameter plane. To take into account the fact 
that at times about 1 — 2 units of Ra the transverse expansion should lead to fast cooling of the 
hot QCD matter we also impose the condition L < Lmax- We performed the computations 
for Lmax = 8 and 10 fm. The difference between these two versions is small. 

4. In Fig. 1 the theoretical Raa obtained for = 0.7, 0.6, and 0.5 for the chemically equilib- 
rium and purely gluonic plasmas is compared to the PHENIX data^^^ on vr'^ production in the 
0-5% central Au + Au collisions at = 200 GeV. The results are presented for radiative energy 





a;'=o.5 





10 20 
Pt [GeV] 



10 



20 



Figure 1: The factor Raa for vr" productfon in the 0-5% central Au+Au collisions at ^/s = 200 GeV for — 0.7, 
0.6, and 0.5. The upper panels are for the chemically equilibrium plasma, and the lower ones for purely gluonic 
plasma. Black line: the total radiative part (quarks plus gluons); red line: the radiative quark part; green line: 
the radiative gluon part; blue line: the radiative (quarks and gluons) plus coUisional, and plus energy gain due to 
gluon absorption. The theoretical curves obtained for L max = 8 fm. The experimental points are the PHENIX 

datallSl. 

loss and with inclusion of collisional energy loss and radiative energy gain. The effect of the 
radiative energy gain on Raa is practically negligible and can be safely neglected. The growth 
of Raa for gluons in Fig. 1 is due to the q ^ g transition which is usually neglected. However, 
it does not affect strongly the total Raa since for ^/s = 200 GeV the gluon contribution to the 
hard cross section is small at pT > 15 GeV. In Fig. 2 we compare our results for a(''' = 0.7, 0.5, 
and 0.4 with the ALICE data^^l for charged hadrons in Pb + Pb collisions at ^/s = 2.76 TeV. 

As can be seen from Figs. 1, 2, the collisional energy loss suppresses Raa only by about 
15-25%. For the equilibrium plasma the data for ^/s = 200 GeV can be described with ~ 
0.6 0.7. The data for = 2.76 TeV agree better with ^ 0.4 -+ 0.5. It provides evidence 
for the thermal suppression of at LHC due to higher temperature of the QGP. 



5. In summary, we have analyzed the data on Raa obtained in the PHENIX experiment on 
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Figure 2: The same as in Fig. 1 for the charged hadrons in Ph + Ph collisions at \/a = 2.76 TeV for = 0.7, 
0.5 and 0.4. The experimental points are the ALICE data^, as in^the boxes contain the systematic errors. 

Au + Au collisions at = 200 at RHIC and in the ALICE experiment on Ph + Ph 

collisions at -^/i = 2.76 TevEl at LHC. Our results show that slow variation of Paa from RHIC 
to LHC supports that the QCD coupling constant becomes smaller in the hotter QGP at LHC. 
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